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Abstract 

In this paper, we propose a new self-triggered consensus algorithm in networks of multi-agents. 
Different from existing works, which are based on the observation of states, here, each agent determines 
its next update time based on its coupling structure. Both centralized and distributed approaches of the 
algorithms have been discussed. By transforming the algorithm to a proper discrete-time systems without 
self delays, we established a new analysis framework to prove the convergence of the algorithm. Then 
we extended the algorithm to networks with switching topologies, especially stochastically switching 
topologies. Compared to existing works, our algorithm is easier to understand and implement. It explicitly 
provides positive lower and upper bounds for the update time interval of each agent based on its coupling 
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structure, which can also be independently adjusted by each agent according to its own situation. Our 
work reveals that the event/self triggered algorithms are essentially discrete and more suitable to a 
discrete analysis framework. Numerical simulations are also provided to illustrate the theoretical results. 


Key words: self-triggered, distributed, consensus, switching topology. 


I. Introduction 

Distributed control of networked cooperative multiagent systems has received much attention in 
recent years due to the rapid advances in computing and communication technologies. Examples 
include agreement or consensus problems [ffl|-[[3]|, in which a group of agents seek to agree upon 
certain quantities of interest, formation control of robots and vehicles @j, 0, and distributed 
estimation m, a , etc. 

In these applications, an important aspect is to determine the controller actuation schemes. 
Although continuous feedback of states have always been used in early implementations of 
distributed control, it is not suitable for agents equipped with embedded microprocessors that 
have very limited resources to transmit and gather information. To overcome this difficulty, 
event-triggered control scheme lf8l- lfT6l was proposed to reduce the controller updates. In fact, 
event-driven strategies for multi-agent systems can be viewed as a linearization and discretization 
process, which was considered and investigated in early papers Ifni . Ifl8i . For example, in the 
paper iflTft . following algorithm was investigated 

m 

x\t + 1) = f(x l (t)) + Ci aij(f (x 3 (t))) (1) 

3 = 1 

which can be considered as nonlinear consensus algorithm. In particular, let f(x(t)) = x(t) and 
Ci = (f fc+1 — t l k ). Then © becomes 

m 

x\t l k+1 ) = x\t\) + (4 +1 - t' k ) (2) 

3 = 1 

which is some variant of the event triggering (distributed, self triggered) model for consensus 
problem. In centralized control, the bound for (t l k+1 — t\) = {t k +1 — t k ) to reach synchronization 
was given in the paper ifTTIl when the coupling graph is indirected and in |fl8l for the directed 
coupling graph. The key point in event-triggered algorithm is the design of a decision maker that 
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determines when the next actuator update should occur. The existing event-triggered algorithms 
are all based on the observation of states. For example, in a typical event-triggered algorithm 
proposed in (S), an update is triggered when a certain error of the states becomes large enough 
with respect to the norm of the state, which requires a continuous observation of the states. 
In addition, self-triggered control Ifl9ll - ll23ll has been proposed as a natural extension of event- 
triggered control, in which each agent predicts its next update time based on discontinuous state 
observation to further reduce resource usage for the control systems. 

Both the event-triggered and self-triggered control algorithms that have been proposed till 
now have some drawbacks in common. First, the resulting system in event-triggered control is 
one that with system delays, especially with self delays, which generally is difficult to handle. 
And the existing analysis for such algorithms are always based on some quadratic Lyapunov 
function, which has very strict restrictions on the network structure. For example, to the best 
of our knowledge, the latest analysis is still restricted to static networks. Second, since these 
algorithms are based on the observation of states, which are generally much complicated and 
untraceable, making it difficult to predict and exclude some unexpected possibilities such as the 
occurrence of zero inter-execution times. 

To overcome such difficulties, we proposed a new ki nd of self-triggered consensus algorithms 
that are structure-based. That is, each agent predicts its next update time based on its coupling 
structure with its neighbours instead of the observation of their states. This can bring several 
advantages. First, since the coupling structure is relatively simpler and more traceable than the 
states, it is relatively easier to handle. Actually, in our algorithm, the occurrence of zero inter¬ 
execution time has been excluded by directly providing a positive lower bound on the update 
intervals. 

Second, although the resulting system of the self-triggered algorithm is one with self delays, 
by some proper transformation, we showed that the system can corresponding to a discrete-time 
system without self-delays. We have proposed the algorithm in both centralized and distributed 
approaches. In the centralized approach, the system can be directly related to a discrete-time 
system without delays, which can be seen as a discrete version of the nominal system. However, 
in the distributed approach, the situation is much more complicated and the system can not 
be directly transformed into its discretized version. However, by some proper indirect transfor¬ 
mation, we showed that the convergence of the nominal system can be reduced to that of a 
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discrete-time consensus algorithm with off-diagonal delays but without self delays. Thus, other 
than using a quadratic Lyapunov function, the convergence of the original algorithm can be 
solved by the analysis of a discrete-time consensus algorithm. This brings another possibility 
that is also considered as a big advantage and a major contribution of this work. That is, we can 
analyze networks with switching topologies, especially stochastically switching topologies. 

The rest of the paper is organized as follows: Section HD provides some preliminaries on 
matrix and graph theory that will be used in the main text. Section Hill provides the self-triggered 
consensus algorithm in both centralized and distributed approach with convergence analysis. 
Section [IV] provide an example with numerical simulation to illustrate the theoretical results 
obtained in the previous section. The paper is concluded in Section [V] 

II. Preliminaries 

In this section, we present some definitions and results in matrix and graph theories and 
probability theories that will be used later. 

For a matrix A = [a r; ], a i? represents the entry of A on the /1h row and yth column, which is 
sometimes also denoted as [A],j. A matrix A = [a tJ ] is called a nonnegative matrix if a l:j > 0 
for all i, j. And A is called a stochastic matrix if A is square and = 1 for each i. Given 

a nonnegative matrix A and 8 > 0, the 8-matrix of A is a matrix that has nonzero entries only 
in the place that A has entries equal to or greater than 8. For two matrix A, B of the same 
dimension, we write A> B ii A — is a nonnegative matrix. Throughout this paper, we use 
^ =1 Ai = A k A k _i ■ ■ ■ Ai to denote the left product of matrices. 

A directed graph Q is defined by its vertex set V{Q) = {iq, • • • : v n } and edge set £{G) C 
{(vi,Vj) : Vi,Vj E V(G)}, where an edge is an ordered pair of vertices in ViQ). If (v t . Vj) e 
£(G), then u, is called the (in-)neighbour of Vj. Mi = {j : ( Vj,Vi ) e £{G)} is the set of 

neighbours of agent i. A directed path in Q is an ordered sequence of vertices Vi, ■ ■ ■ , v s such 
that (v i: v i+ i) e £{Q) for % = 1, • • • ,s — 1. A directed tree is a directed graph where every 
vertex has exactly one in-neighbour except one ,called the root , without any in-neighbour. And 
there exists a directed path from the root to any other vertex of the graph. A subgraph of Q is 
a directed graph Q s satisfying V(Gs) C V(f?), and £{Q) C £{G). A spanning subgraph of Q is 
a subgraph of Q that has the same vertex set with Q. We say Q has a spanning tree if Q has a 
spanning subgraph that is a tree. 
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A weighted directed graph is a directed graph where each edge is equipped with a weight. 
Thus, a graph Q of n vertices corresponds to an n x n nonnegative matrix A = a^], called 
the weight matrix in such way that ( Vj,Vi ) £ £(G) if and only if a tJ > 0. On the other hand, 
given an n x n nonnegative matrix A, it corresponds to a weighted directed graph G(A) such 
that G(A) has A as its weight matrix. Using this correspondence, we can introduce the concept 
of 5-graph. The 5-graph of G is a weighted directed graph that has the 5-matrix of G' s weight 
matrix as its weight matrix. A graph G has a 5 spanning tree if its 5-graph has a spanning tree. 

From the weight matrix A = [a^] of a graph G, we define its graph Laplacian L = [l t j] as 
follows: 


hj — 


'Yhk^i 0j iki — -i > 


-aij , i j. 


Thus L has zero row sum with l ri > 0 while < 0 for i,j = 1,2, ■ • • ,n and i j. And 
the nonnegative linear combination of several graph Laplacians is also a graph Laplacian of 
some graph. There is a one-to-one correspondence between a graph and its weight matrix or its 
Laplacian matrix. In the following, for the sake of simplicity in presentation, sometimes we don’t 
explicitly distinguish a graph from its weight matrix or Laplacian matrix, i.e., when we say a 
matrix has some property that is usually associated with a graph, what we mean is that property 
is held by the graph corresponding to this matrix. For example, when we say a nonnegative 
matrix A has a spanning tree, what we mean is that the graph of A has a spanning tree. And 
it is of similar meaning when we say that a graph has some property that is usually associated 
with a matrix. 

Let {U, T ', P} be a probability space, where Q is the sample space, T is a-algebra on 0, and 
P is the probability on T. We use E{-} to denote the mathematical expectation and E{■ | J 7 } the 
conditional expectation with respect to J 7 , i.e., 1E{-I^T 7 } is a random variable that is measurable 
with respect to T. 

Definition 1 (adapted process/sequence): Let {A k } be a stochastic process defined on the 
basic probability space {U, J 7 , P}, and let {J 7 *,} be a filtration, i.e., a sequence of nondecreasing 
sub-a-algebras of T. If A k is measurable with respect to J 7 /,, then the sequence {A k . 5F k ) is 
called an adapted process/sequence. 

Remark 1: For example, let X (cu) be a random variable defined on O, then {E{X(cu) \J r k }. J r k \ 
is an adapted sequence. 
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III. Self-Triggered Consensus Algorithm 


In a network of n agents, with x r G E being the state of agent i, we assume that each agent’s 
dynamics obey a single integrator model 

£i(t) =Ui(t), i — 1,2, ■ ■ ■ ,n, (3) 

where Ui(t) is the control law for agent i at time t. In consensus algorithm, the control law is 
usually given by( flU, |E3, 11301 ) 

n 

M 1 ) = -'52 kjMt) - x i(t )] = -52 ( 4 ) 

je Mi j =i 

Thus the consensus algorithm can be written as 

n 

Xj(t) = ~y ^Jij[xj(t) - Xj(t )], (5) 

3 = 1 

or in matrix form as 

x(t) = —Lx(t), (6) 

where x(t) = \x\ (t), ■ ■ ■ ,x n (t)] T G M n is the stack vector of the agents’ states. In event/self- 
triggered control, the control law u(t) is piecewise constant, i.e., 

u(t) = u(t k ), t G [t k ,t k+ 1 ), (7) 

where {tk} is the time sequence that an update of u (•) occurs. The central point of such algorithm 
is the choice of appropriate {tk} such that some desired properties of the algorithm such as 
stability and convergence can be preserved. This can be done in a centralized approach or a 
distributed approach, both of which will be discussed in the following. 

A. Centralized Approach 

In this subsection, we first consider centralized self-triggered consensus algorithm. The se¬ 
quence of the update time is denoted by: to, t\, t 2 , ■ ■ ■, then the self-triggered algorithm has the 
form: 

x{t) = —Lx(t k ), t G [t k ,t k+ 1 ). (8) 

Denote A t k = t k+1 — t k for k — 0,1, • • •. We have the following Theorem. 
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Theorem 1: If the graph G(L) has a spanning tree and there exists 5 G (0,1/2) such that 
S /I m ax < Af fe < (1 — S)/l max for each k, where Z max = maXj{4}, then the algorithm will 
achieve consensus asymptotically. 

Proof: From the algorithm, for each t G (4,4+ 1 ], since 

x(t) = -Lx(t k ), 


we have 

x(t) = x{t k ) -{t- t k )Lx(t k ) = [I - (t - 4)L]x(4). 
Particularly, let t = 4 + 1 , we have 


x(4+i) = (/ - At k L)x(t k ) = A k x(t k ), 


where = I — AIt is easy to verify that A is a stochastic matrix for each k. In fact, 
since 5 / /max ^ — (1 ^)/^max? 

1 ^4^max ”7 /max 

^max 

For i ^ j. 


[A k \ij — —At k lij > 0. 


It implies that [A k ]ij is positive if and only if —l l:] is positive, and the positive elements of [A k \ 
is uniformly lower bounded by a positive scalar since both At k and — is uniformly lower 
bounded by a positive scalar. Furthermore, Y^j=i[A-k\ij = 1 — At k , hj = 1. By a standard 
argument from the theory of products of stochastic matrices, there exists i*gR such that 


lim Xi{t k ) = x*. 

k— >+oo 

On the other hand, it is not difficult to verify that the function maxi{a+(t)} is nonincreasing and 
minj{xi(f)} is nondecreasing. Thus, 

lim xAt) = x*. 

£->-+00 


Using a similar argument as that in Theorem Q] it is not difficult to propose a centralized event- 
triggered algorithm on networks with stochastic switching topologies. For example, consider the 
following consensus algorithm: 

x(t) = —L k x(t k ), t G [4,4+i), (9) 
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where L k = [/■'■] is uniformly bounded in k and l k YAX = maXj{/ ? fc J e [/ m i n , /max] for some 
0 < /min < /max- Before provide the next theorem, we first summarize the following assumption. 
Assumption 1: (i) There exists 5 6 (0,1/2) such that At fc 6 [5//^, (1 — S)/l^ J; 

(ii) {L fc , is an adapted random sequence, and there exists h > 0, 8' > 0 such that the 
graph corresponding to the conditional expectation ^{Y ^= k +1 L m \ Fk} has a (^'-spanning 
tree; 

Now we have 

Theorem 2: Under assumptiorQ] the self-triggered consensus algorithm © will reach consen¬ 
sus almost surely. 

Proof: Similar as in the proof of Theorem Q] we have x(t k+ 1 ) = A k x{t k ) with A k = 
I — At k L k . Since {A t k } is a deterministic sequence, {A kl 2F k } is also an adapted sequence, and 

k+h k+h 

E{ A ™\Fk} = hi- E{ ^ At m L m \F k } 

m=k+l m=k+ 1 

r- k + h 

> h5I-f- E{ ^ L™\F k }, 

^max , . , 

m=k+l 

where [L™] i:j = [L 171 ]^ for i f j and [L™] u = 0 for all i. Thus, ^{Y^m= k +i A rn\F k } has a 
^"-spanning tree with 8" = 8'8/l max > 0. From Theorem 3.1 of |[28l1 . the sequence {x(i fc )} will 
reach consensus almost surely, thus x(t) will also reach consensus almost surely. ■ 

As corollaries in special cases, it can be shown that almost sure consensus still holds if we 
replace item (ii) in Assumption Q] with one of the following more special ones. 

(ii’) {L k } is an independent and identically distributed sequence, and there exists 8' > 0 such 
that the graph corresponding to the expectation E L k has a (^'-spanning tree; 

(ii”) {L k } is a homogenous Markov chain with a stationary distribution n, and there exists 
5' > 0 such that the graph corresponding to the expectation with respect to it, L k , has 
a ^'-spanning tree; 

B. Distributed Approach 

In this section, we discuss the distributed event-triggered consensus algorithm. 

Let {t k } k = 0 be the time sequence such that t k is the /. th time that agent i updates its control 
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(10) 


law. Then the distributed self-triggered algorithm has the following form: 

n 

xft) = - y2 = 1,-*- ,77,, 

3 =1 

where kj(t) = max{fc : t 1 * < f}. 

Denote Atf = t\ +1 — t\. Then we have the following Theorem. 

Theorem 3: If the graph G(L) has a spanning tree and there exists c), e (0,1/2) for i = 
1,2, •• • , n such that AG [(//)„;, (1 — Si)/la], then the distributed event triggered algorithm 
can reach a consensus as t —> oo. 

The proof will be divided into two steps corresponding to two lemmas. Now we give and 
prove the first lemma. 

Lemma 1: If there exists i*eR such that 


lim x{t\) = x*, i — 1, 2, ■ ■ ■ , n, 

k —^oo 

then the distributed event-triggered algorithm will reach a consensus, i.e., 

lim Xi(t ) = x*, i — 1, 2, ■ ■ ■ , n. 

t—>oo 


Proof: By the definition of the algorithm, for agent /', 

n n 

Xt{t) = - (<) ) = - li i\- x 3 (tf® ) - Xi ■ 

3 = 1 

Since lim^oo kj(t ) = +oo, 


Thus, 


lim xAt) 

t—> OO 


n 

Y. <«(*•-*•) = o. 


lim | xft) — x* 

t—yoo 


< lim \xi(t) - xft^f + lim | xft^) - x*\ 

t—> OO t—>oo 

< lim Atf^ max liAs)! + lim IxAt^^) — x* 

t^oc selt^i] t -*’° 

< , lim max |xj(s)| + lim I xAt k f^) — x* 

la t] 


Now we are to provide and prove the second lemma. 
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Lemma 2: Under the assumption in Theorem 0 there exists x* such that 

lim xft 1 *) = x*. 

k —^OO 

Proof: Let {t k } be the time sequence that the /cth update occurred in the network, i.e., 
{t k } = U” =1 {ff }, and A t k = t k+1 — t k . In case for that t k = tf = tf for some some k, k', 
k" and i f j, we consider the update of agent i and j at time t k as one update of the network. 
Now, construct an auxiliary sequence {y(k)}, where y(k) = [yi(k),y 2 (k), ■ ■ ■ ,y n (k)] T G M" 
with yfk) = Xi(t^ tk '), i.e., yfk) is the latest state that agent i uses in its control law at time 
t k . Now consider the evolution of {y(k)}. 

Case 1. Agent i does not update at time t k+ i. In this case, by definition, 

yfk + 1) = yi(k)- 

Case 2. Agent i does update at time t k+1 : In this case, tf‘ ( ' tk+1> = t k+l , assume tf' < ' tk ' > = t k _ d . k be 
the last update of agent i before t k + 1 , with d ik > 0 being the number of updates occur at all 
other agents Xj, j i, between the two successive updates of agent i, i.e., in the interval 
By definition yfk) = • • • = yfk - d ik ), and xft) = - YTj =1 hjVjik') for 
t G (t k ',t k '+i) for all k'. Since [t^ tk \t^ th+l) ) = U^ =fe _ d .Jf TO , t m +i). 


10 


Ui(k + 1) = Xi(t k+1 ) = Xiit'l^) + / 

J t k i^k) 

djk pk—dik+m+l 

= x^dj + £ / 

771=0 tk ~ d ik+ rn 
d*ik Tl 

%i{th—dik') ^ ^ dik+m ^ ^ ^ijVj ^zfc “1“ 

m=0 j=l 

dik 

2/z(^ ^zfc) ^ ^ ^tk—dik~\-mliiyi(k 6^z/c “1“ 
ra=0 

dik Tl 

^ ^ At k _d ik -\- m ^ ^ hjVj(k d^ k -|- Tfl) 


m 


m=0 j=l,j& 

d{k dik Tl 

2/i(fc) ^ ^ ^tk— ^fc+m^zz?/z(^) ^ ^ ^tk—dik+m ^ ^ &ik H - 


m) 


m =0 
dik 


m =0 

dik Tl 


j=l,j^i 


(1 ^zz ^ ^ ^^k—dik~\-m)yi (^) ^ ^ ^ ^ ^^k—dik~\-Tn^ijyj ^ik “1“ ^7») 

m=0 m=0 j=l,j^i 

dik Tl 

= (1 - A t^ tk) lii)yi(k) - ^2 'Yl At k-d ik +mhjyj(k - d ik + m) 

m=Oj=l,jjti 

dik n 

171=0 j = 1 

where ct° (fc) = 1 - A t\ i{ - tk) la, a%(k) = 0 for m = 1, 2, ■ ■ ■ , d ik and a%(k) = -At k _ m l iy 
It is easy to verify that a% (k) = 1 — A t^^la > 1 — (1 — <5*) = St, a™j = —A t k _ m kj > 0, 
and 

dik n 

EE<w = L 

772=0 j = 1 


Besides, by the assumption 8i/lu < At\ < (1 — 8i)/la for each i, k, it is clear that 
in the interval [Z*,Z* +1 ], only finite updates occur for agents x 3 with j ^ i, and the 
number of updates is uniformly upper bounded, i.e., there exists r > 0 independent of 
i, k, such that d ik < t. For example, let <5 min = liiin,; {(),//„}, and 5 max = max, {(),//„}, 
then 5 min < A t\ < 5 max for all i, k, and on each time interval of length m<5 min , at most 
m + 1 updates occur for each agent. Let h' be the smallest positive integer satisfying 


11 


h'5 min > 5 max , then on each time interval of length 5 max , at most h' + 1 updates occur for 
each agent. Thus on the interval (t^,t^ +1 ), at most (n — 1 )(/T + 1) updates occur for all 
agents j with i ^ j. Pick r = (n — 1 ){h' + 1), then < r for all i, k. 

If agent i does not update at time tk+i, i.e., case 1, we define a™(k) = 0 for all m = 0,1, • • • , r 
and i, j — 1 , 2, • • • , n except for a° (k) = 1 , 

If agent i updates at time tk+i, i.e., case 2, we define a^(k) = 0 for all = and 

m = dik + 1, • • • , r when < r. 

For both cases, we can give a uniform iterative formula for y(k) as: 

r n 

y(k+l) = J2^2 al ij(k)yj(k-l), ( 11 ) 

i=o j =i 

where > 0 for all i,j , L k, and 


EE 

1=0 j =1 


%-(k) = 1, 


for all i, k. Denote 5 = min,{5j} > 0, we have «■((/;:) > 5 for all i, k. 
Construct a new matrix B(k) = [bij(k)\, where 


t>u(k) = al{k ) 


( 12 ) 


and 

r 

b ij( k ) = i f 3 (13) 

1=0 

Obviously, B(k) > 51, and b l3 (k) = A> 5 m i n kj if agent i updates at time tk+i- 
If there exists h > 0 such that on each interval [f fe ,f fc+fe ] all agents update at least once, then 
Y^m=k B{m) > —5 min L for each k. Since <3(L) has a spanning tree, there exists 5' > 0 such that 
Y?m= k B (m) has a ^'-spanning tree for each k. From Corollary [2] in the appendix, the sequence 
{y(k)} will reach a consensus, i.e., there exists i*6l such that 


lim yt(k) = x*. 

k —^oo 

This is equivalent to 

lim £j(tf) = x*. 

k—yoo 


At last, we give some hint on how to calculate the quantity h mentioned above. It is easy to see 
that each agent updates at least once on each time interval of length () rnax , since A < c) max for 
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all i, k. And in each time interval of length 3 max , there are at most n{h! + 1) updates of all the 
agents, where h! is defined as before. Pick h = n ( h' + 1), then for each k, each agent updates 
at least once in the time interval [t k ,t k+h \. ■ 

Remark 2: It is not difficult to see that this algorithm is applicable to leader-follower networks, 
in which the leader receives no information from others. For example, if agent i is the leader, 
then la = 0 and the update time interval A= +oo by definition, which mean the leader i does 
not need to update its state. Generally, if agent i wants a long update time interval, all it needs 
to do is to choose small coupling weights so that is small, and this can be done independent 
of all other agents. 

From the above analysis, it is not difficult to see that the distributed event-triggered algorithm 
can also be applied to networks with switching topologies. For example, at each update time 
agent i may reset its coupling weights Uj such that the algorithm can be written in the form: 

n 

Xi(t) = - l ij x j ( tk / {t) )> t e [t k , f • +1 ), i = 1, 2, • ■ ■ , n, (14) 

3 = 1 

where <5*//^ < A t\ < (1 — Si)/Ik. One of the great advantage of this algorithm is that each 
agent can adjust its update time interval independently based on its own need. If agent i wants 
a longer update time interval, then it may decrease the coupling weight, otherwise, increase the 
coupling weight. For example, each agent i may store a scaling factor ef such that 

J,‘ = £,%• (15) 

By a similar analysis as in the case of fixed topology, we can prove 

Theorem 4: Suppose that e min and e max are two positive number satisfying 0 < e min < e max 
and for all i, k, e min < e* < e max . If the graph Q{L) has a spanning tree, then the distributed 
event-triggered consensus algorithm with switching topology (fl4l) and weight updating rule (IT5l) 
can reach consensus as t —» oo. 

Finally, we investigate event-triggered consensus algorithm in networks with stochastically 
switching topologies in which each agent i updates its coupling weights independently at the 
same time it updates its state. 

Given 0 < a < b, and let Si(a, b) = {s = [si, ■■-,«„]: s* e [a, b ], .Sj < 0, j ^ i, YTj=i s j = 

0}. The event-triggered algorithm can be formulated as follows: 

n 

x i(t) = ~ ~ lk ij x i (*? W )> i 6 ( 16 ) 

3 = 1 
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where L k = [l k { , • • • ,l k n ] G S' ? (a ): . /;,) for some given 0 < a, < b t , and {is an independent 
and identically distributed sequence. Before stating the convergence result, we summarize the 
following assumption. 

Assumption 2: (i) There exist Si G (0,1/2) such that At* G [Si/l\ (1 — <5j)//|)]; 

(ii) There exists <5 > 0 such that the graph with Laplacian being E L k has a ^-spanning tree, 
where L k = [T^ T , • • • , T^ T ] T is a matrix with its zth row being L k \ 

Theorem 5: Under Assumption [2} the event-triggered consensus algorithm (fl6l) will reach a 
consensus almost surely. 

Proof: The algorithm (fl6l) can be reformulated as: 

n 

Xi(t) = -£fe(^), t e [4,4+1), (17) 

3 = 1 

where /•) = Let L k = [Z4], then generally the sequence {L k } is not an independent 

sequence since for each k, l k - = /^ +1 if no update of agent i occurs at time t k+1 . However, 
similar to the analysis given above, it can be seen that L k is independent of L k ' for k' > k + N, 
where N = n(h'+ 1) and h! is the smallest positive integer satisfying h' min* {Si/hi} > max i {(l — 
Si)/ai}. 

Similarly, we have 

N n 

y(k+l)=^2^2° l ij (k)y j (k-l), (18) 

1=0 j =1 

where a A (k) > 0 for all i,j, l, k, and 

N n 

1=0 j =1 

for all i, k. Since a^-(fc) = —A t k _il k ~\ we can see that A°(k), ■ ■ ■ ,A N (k) are independent of 
A°(k'), ■ ■ ■ ,A N (k'), when k' > k + 2 N. On the other hand, in each time interval ( t k ,tk+N ), 
each agent there at least updates once. Thus we have 

k+N 

E B(k)>5'L 0 , 

m=k-\- 1 

where B(k) is defined as before and S' = min, ; fc {Af/} = min v {} > 0, and [Z 0 ]^ = 
— [E L k ]ij for i f j and [L 0 \a = 0 for all i. This implies that Y^=k+\ E B(k) has a ^"-spanning 
tree with 5" = 55' > 0. The conclusion follows from Corollary [3] in the Appendix. 
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IV. Numerical Simulation 


In this section, we provide an example to illustrate the theoretical results in the previous sec¬ 
tion. We use the distributed event-triggered algorithm in networks with stochastically switching 
topologies as considered in Theorem [5] 

Consider a network of four agents with 

Lie {[1,-1,0,0],[i,o,o,-i]}, 

L\ = {[- 1 , 1 , 0 , 0 ], [ 0 , 1 ,- 1 , 0 ]}, 

L k 3 e {[0,-1,1,0], [o,o, 1,-1]}, 

Z 4 fc e {[o,-i,0,1], [0,0, 

and each agent selects its coupling weights using a uniform distribution. We choose 5 t = 0.1 
for all the agents. The next update time is randomly chosen from the permissible range. It is 
not difficult to verify that the conditions in Theorem [5] are satisfied. The simulation results are 
provided in Fig. Q] with the initial value of the four agents being randomly chosen. It can be 
seen that the agents actually reached consensus. 

V. Conclusion 

In this paper, we proposed a new structure-based self-triggered consensus algorithm in both 
centralized approach and distributed approach. Different from existing works that used a quadratic 
Lyapunov function as the analysis tool, we reduce the self-triggered consensus algorithm to a 
discrete-time consensus algorithm by some proper transformation, which enables the application 
of the product theory of stochastic matrices to the convergence analysis. Compared to existing 
work, our method provides several advantages. First, each agent does not need to calculate 
the system error to determine its update time. Second, it provides explicit positive lower and 
upper bounds for the update interval of each agent based on its coupling weights, which can 
be adjusted by each agent independent of other agents. We also used our method to investigate 
networks with switching topologies, especially networks with stochastically switching topologies. 
Our work reveals that the event/self-triggered algorithms are essentially discrete and thus more 
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Fig. 1. Self-triggered consensus in networks of four agents with stochastically switching topologies. 


suitable to a discrete analysis framework. And it is predictable that more results can be obtained 
for self-triggered algorithms under our new analysis framework in the future. 

Appendix: Discrete-time consensus analysis 

The appendix is devoted to the discussion of discrete-time consensus algorithms in networks 
of multi-agents with time delays based on the product theory of stochastic matrices and provide 
some results that is used in the main text. 

A. Preliminaries 

This subsection will provide some preliminaries on stochastic matrices that will be used in 
the proof of Theorem [6] in the next subsection. 

A nonnegative matrix A is called a stochastic indecomposable and aperiodic (SIA) if it is a 
stochastic matrix and there exists a column vector v such that lim^oo A k = e n v T , where e n is 
an n-dimensional column vector with all entries 1. A sequence of stochastic matrices {A k } is 
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ergodic if there exists a column vector v such that lini/,^ oc n^Li -'U = e„'<; T . A is called 5-SIA 
if its 5-matrix is SIA. A stochastic matrix A is called scrambling if for any pair (i, j), there 
exists an index k such that both a ik > 0 and a :jk > 0. Similarly, A is called 5-scrambling if its 
(5-matrix is scrambling. For a matrix A, diag(A) refers to the diagonal matrix formed by the 
diagonal elements of A, i.e., [diag(A)]^ = [A]a for all i and [diag(A)]y = 0 for i ^ j. And we 
define diag(A) = A — diag(A). 

We use Pl'lJ 7 } to denote the conditional probability with respect to J 7 on a probability space 
{0, P, T 7 }, which can also be expressed into the conditional expectation of an indicator function, 
i.e., Pj'lJ 7 } = E{1{.}|J 7 }, where 1/.} is an indicator function which is defined as 


!s(z) = 


1, x e S, 

0, x qL S. 

Thus, P{• I^T 7 } is also a random variable that is measurable with respect to T. 

In the following, we will introduce some lemmas that will be used in the proof of the main 
results. 

Let A = [a* j] be an n x n stochastic matrix. Define 


A (A) = max max | a^j — a i2 j |. 

j ii,»2 

It can be seen that A (A) measures how different the rows of A are. A (A) = 0 if and only if 
the rows of A are identical. Define 


A(A) = 1 — min ) min{a lu , a l2j }. 

ii,i2 ^— J 

j 


Remark 3: It can be seen that if A is scrambling, then A < 1, if A is 5-scrambling for some 
5 > 0, then A(A) <1 — 5. 

Lemma 3: lf24l For any stochastic matrices A 1 , A 2 , • • •, A k , k > 0, 

k 

A(A 1 A 2 ---A fc )< J]A(A). 

i=l 

Remark 4: From Lemma (3} it can be seen that for a sequence of stochastic matrices {A k }, if 
there exists a sequence {ki} such that ki < k i+ 1 , lim^oo ki = oo, and fl^i Adlm^Wi An) = 0’ 
then linifc^oo A(f|^ =1 A m ) = 0, i.e., {A/,.} is ergodic. Particularly, from the property of scrabling- 
ness, if there exists 5 > 0 such that there are infinite matrices in the sequence {Jlmi/c.+i An} 
that are 5-scrambling, then {A k } is ergodic. 
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Lemma 4: 112411 . [261 Let A lt A 2 , ■ ■ ■, A k (repetitions permitted) be r x r SIA matrices with 
the property that for any 1 < k\ < k 2 < k, YliL kl A is SIA. If k > n(r), then nf=i A ’ s a 
scrambling matrix. Here n(r) is the number of different types of all r x r stochastic matrices. 
Lemma 5: ll26ll Let A\, • • •, A rn be n x n nonnegative matrices, let 

A\ A 2 ■ ■ ■ A m 

0 ••• 0 


D = 


mnXmn 


let 


0 

0 0 

I 0 
I 0 
0 I 

0 0 

and let Mk — D + Mq for any k e {1, 2, • • • , m — 1}. Then if GiYl'iLi A) contains a spanning 
tree, then Q{M k ) contains a spanning tree with the property that the root vertex of the spanning 
tree has a self-loop in G(M k ). 

r I 0 ••• 0 


M 0 = 


0 0 
0 0 
0 0 
0 0 
I 0 


mnxmn 


Remark 5: Actually, since Mq = M,"' 1 = 


I 0 


0 


for all k > m — 1, thus Lemma 


:>ased on this lemma. 


/ 0 ••• 0 _ 

[5] holds for all k £ N. This is important for our further results 

Remark 6: If the condition is h as a ^-spanning tree for some 5 > 0, then the 

conclusion becomes that G(M k ) has a ci'-spanning tree whose root vertex has a self-loop, where 

5' > 5/m. 

Lemma 6: lf25l Let A be a stochastic matrix. If G(A) contains a spanning tree with the 
property that the root vertex of the spanning tree has a self-loop in G{A), then A is SIA. 

From Lemma [5] and [6] we can have the following Lemma. 
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Lemma 7: Let A\, • • •, A l m , i = 1, • • • . k be n x n nonnegative matrices. Let 


D, 


M A i ... A i 

0 0 ••• 0 


0 0 ••• 0 


mn x mn 

and M 0 being defined as in Lemma [5] Then if G(52i =i zLj - 1 ) has a (^-spanning tree for some 

5 > 0, then the product n*Li(A + M)) has a ^'-spanning tree for some 0 < 5' < 5. And the 
root of the spanning tree has a self-loop, thus fl !'-1 (A + M 0 ) is ASIA. 

Remark 7: It is obvious from Lemma [5] that if there exists e > 0 such that D\ > e(M 0 + A), 
then nti A is ASIA for some 0 < 8' < 8. This is the case that will appear in the proof of 
Theorem [6j 

Proof: First, we have 

k k k 


J](M 0 + A) + + (19) 

i= 1 i =1 i= 1 

where the second inequality is due to the fact that Mg A > Di for any j > 0. And it is not 
difficult to verify that the first block row sum of 1), is preserved in D,M'[f^ for % — 1, • ■ ■ , k. 
Thus the first block row sum of Y^l=\ A^fA* equals that of Yli=i A. he., Yli=i Sj=i A'- This 
means AAA 1 ) has a 5-spanning tree. From Lemma [51 Q(Mq + ^ =] AMq - 1 ) has a 

Aspanning tree for some 0 < 5' < 8, and the root has a self-loop. This is also true for the 
product ntr(M 0 + Di) due to (fl9l) . Thus, Yl^ = i(M 0 + Di) is ASIA from Lemma [6] The proof 
is completed. ■ 

From Lemma [7] we can easily obtain the following corollary. 

Corollary 1: Let A\, ■ ■ •, A l m , i = 1, 2, • • • , be uniformly bounded n x n nonnegative matrices. 
Let Di and M 0 be defined as in Lemma [7] If there exists k such that AX)!=i 1 A)) has a 
(^-spanning tree for some 8 > 0, then for each k' > k, the product Yl^ =1 (Di + M 0 ) is 8'{k')- SIA, 
where 0 < 8'{k') < 8 depends on k'. 

Proof: Since A *■ > 0, if has a Aspanning tree for some 8 > 0, then 

for each k' > k, G(Yli=i X);=i A) also has a (^-spanning tree. The conclusion is obvious from 
Lemma [7] ■ 
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Remark 8: If A\, A l m , i = 1,2,--- , is a random sequence, then it is obvious from 
Corollary Q] that the event has a 5-spanning tree is contained in the event 

nf^A + Mo) is 8'{k')- SIA. 

Lemma 8 (second Borel-Cantelli lemma KTTH): Let {A, n > 0} be a filtration with A = 
{0, f2} and |X n , n > 1} a sequence of events with X n e J 7 ,,. Then 

+oo 

{ X n occurs infinitely often } = | ^^P{X n |A-i} = +ooj, 

n= 1 

with a probability 1, where “infinitely often” means that an infinite number of events from 
{X n }n=i occur. 

Lemma 9: Let {Ah, Xu) be an adapted random sequence. If for a sequence {Xf .} defined as 
Xk = {Ak G Sk} for some given set Sk there exists 8 > 0 such that the conditional probability 
P{X m+ i|-A>} > 8, then for each m, h > 0, 

P{X m+1 , • • • , X m+h \X m } > 8 h . 


Proof: By the definition of conditional probability, we have 





= 

E|l.V m+ i ' ' 

1 J 

= 

E{lx m+1 • • 

‘ X rn+h \J~m+h- 

> 


1 ^~LYl J 

> 

<5 2 E{lx m+1 

^ Xm+ h -2 1 } 


> 


> 8 h . 


}|AA 


Lemma 10: (Lemma 5.4 in 11281 ) Let |f2, X, P} be a probability space, and / be a random 
variable with 0 < / < 1. If for a a-algebra J' C 1, a set S G T’ with PIS'} > 0, IE{/1J 27 ’ 7 } > 8 
holds on S for some 8 > 0, then we have 

P{/ > t\F) > 


20 


on S and particularly, 
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on S. 


B. Discrete-time delayed consensus analysis 

In this section, we will discuss discrete-time consensus algorithms in networks with switching 
topologies and time delays. We first consider general stochastic processes described by an adapted 
sequence and establish a most general result. Then we use it to obtain two corollaries under two 
special cases that will appear in the main text. 

Consider the following discrete time dynamical systems with switching topologies: 

r n 

Xi(k+ 1) = 'Y^^2,a\ j (k)x j (k- l), i = (20) 

1=0 j =1 

where k is the time index, x(k) = [x\(k), ■ ■ ■ ,x n (k)] T € R" is the state variable of the system 
at time k, r > 0 is the bound of the time delay, where r = 0 corresponds to the case of no time 
delay. In the following, we always assume that a° d (k) > 0, a l u (k) = 0 for l ^ 0, o.C(fc) > 0 for 
each i, j, l, k, and 

r n 

= 1 

1=0 j =1 

holds for each k and i. 

Define A l (k) = [a l i:j (k)\, and B(k) = [bij(k)] = YZ=o-A!(k). 

We have the following Theorem. 

Theorem 6: Let {A°(k), ■ ■ ■ , A T (k), lF k } be an adapted process. If there exists <5 > 0, h > 0 
such that B(k) > 51 and the conditional expectation iB{k)\5F m } has a 5-spanning 

tree, then the system (l 20 l) will reach consensus almost surely. 

Remark 9: When r = 0, the system has no time delays,and Theorem [6] reduces to Theorem 
3.1 of lf28ll . From this point of view, this result can be seen as an extension of that in lf28l . 

The proof of this theorem will be divided into several steps. 

First, we will prove the following lemma, which shows that the conditional expectation of 
a spanning tree implies a positive conditional probability of existence of a spanning tree for a 
longer length of the summation. 


21 


Lemma 11: Let {A k , lF k } be an adapted random sequence of n x n stochastic matrices, if 
there exists 8 > 0 such that for each m, E{A m+1 \T rn ) has a 8 -spanning tree, then there exist 
h > 0, 0 < 8' < 8 such that 

m-\-h 

P{ ^ A k has a ^'-spanning tree \J r m } > 8'. (21) 

k=m -\-1 

Proof: Denote Sp(n) the number of different types of spanning trees composed of n vertices. 
For some fixed m, we decompose the probability space Q with respect to T m , i.e., Sm ’ G T m , 
fl = US 1 ™^, i\ — 1, 2, ■ ■ ■ , si where s± < Sp(n), and Sm fl Sm = 0 for i f j such that on 
each Sm , E{ A m+1 | T m } has a (fixed) specified spanning tree. Denote E { A rn+ 1 \ IF, A (,,> the 
restriction of E {A m \ T m } on Sm\ i.e., 


^4.77i-|-l | J~m\ o(*i) E{.A m _|_i | Trn} X ljn). 

o?n 

Let ST, 5 (E{A m+ i | 8F m } An)) be a 5-spanning tree contained in E{A m+ i | Tm) ah), which 
is arbitrarily selected when more than one choices are available. For each Sm\ we pick an 
edge (j,i) e £(ST 5 (E{A m+1 | A m } s (n))), and let S^ = {[A m+1 ]^ > 8/2} n St\ Then 
S ( ///! G F m+ 1 , and from Lemma ITOl 

r m ] > | 

holds on Sm ^ ■ 

Similarly, we decompose each S%/h with respect to lF m + i, be., S ^, = U^S^ 1 / 2 / with S^/// G 
T m +i and fl for z 2 ,i' 2 = 1, 2, • • • , s 2 , *2 7 ^ T 2 , where s 2 < Sp(n) depends on the 

index i\ such that E{A m+2 I lF m+ \} An^) has a specified (5-spanning tree on each S^+i • 

For each S %$, we pick an edge (j,i) G £(ST, 5 (E{A m+2 | J m+ i} (il , 2 ))) and let 5'^ 2 ) = 

°m+l 

{[A m+ 2 ]i J > 8/2} n then S^ G Jbn +2 and from Lemma OH 

p{s£$1 r m+1 } > s - 

holds on S%/1\\ 

Continuing this process, we can get a sequence: 


c(il) qi.h) q(ili 2 ) m # # 

5 1 5 ^m+2 j * * * ? 1 J 

such that for Z = 1,2,--- , k, S G T m+i , E{A m+/ | .F m+ i_i 

has a (fixed) spanning tree, and 


} 5 ( iiG- 

m+t —1 


sS?'" i,) = {[A MI ]«>«/2}ns”^, 


( 22 ) 
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where the edge (j,i) E £:(ST 5 (E{yl m+i | Fm+i-i} ^-h)))- 

°m+l -1 

For any fixed sequence i u i 2 , ■■■ ,i k , C S^X\ k) C ^ ^ S mi i- We 

choose the edge (j, i) in ( 1221 ) for each in such way that if ST, 5 (E{ 2 l m+ i | F m +i-\} ainh -n)) 

b m+l-l 

hasn’t appeared in the sequence, ST, 5 (E{A m+ i | .F m } (il )), ■ • •, ST 5 (E{A m+ i_i | Tm+ 1 - 2 } Anh-n-i)), 

b m+l-2 

then we choose (j,i) E £(ST«5(E{2l m+ i | .F m+ ;_i} ^i^-n))) arbitrarily. Otherwise, we choose 

b m+l-l 

(. j,i ) E £(ST,5(E{2l m+ i | J r Tn+ {_i}~(i 1 i 2 ...ij))) which hasn’t been chosen before when 

b m+l -1 

ST^EjAn+i | F m+ i-i} -(iiia-i,)) appears in the sequence 

°m+l-1 

ST 5 (E{A m+ i I JSn}^!)), ■ ■ - , ST 5 (E{A m+ i_i I T m+ i- 2 } a (hi 2 -ii-i)) if there still exists such an 

m b m+l-2 

edge. 

Since there at most Sp(n) different types of spanning trees, and each spanning tree has n— 1 
edges, it is obvious when k — (n — 1) Sp(ra), the sets 

k 

St+k" lk> ^ {J 2 Am +i has ^/2-spanning tree}. 

1=1 

This is because for each fixed sequence i\, ■■■, i k , at least one type of spanning trees has 
appeared at least n — 1 times, thus by the edge chosen strategy, each of its edges has been 
chosen at least once. 
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So we have 


iE Am+i has 5/2-spanning tree | F m } 


1=1 


TP/1 I I 77 \ 

— ^ l u ii,i2,- 3k°m+k I r ™f 


E{ Y, 1 


*1,*2,'" i*fc 


<,(*1*2—^) 

°m+k 


F rn } 


E{E{ V 1 f i i i 2' i k) I ^m+fc-l} | F rn} 
Z - J ^m+fc 




— O ^ v 1 q(*1*2*"*A!) I 3 " 

Z . Z ' . Vfc-1 


*1,12,— >*fc 


0 

Z . . - °m+fc-l 

*1 5*2)’” ,*fc-l 


m} 


> 


> i 7 ”} 

^ Z -^ °m|l 


*1 




*1 


V ’ 

The inequality from (1241) to (l25l) is due to the fact 


(23) 


(24) 

(25) 

(26) 


P{^+ 2 fc " 4) I Frn+k- 1} > ^2 
on E+fc-i fc) e Frn+k—h which implies 

5 

E{1 o( i l i 2 "" i fc) | F rre+fc—1/ > 7rlc(*i*2—*fc)> 

°m-\-k A °m-\-k — 1 

The equality from (1251) to (1261) is due to the fact that S'Efc-i'” ' = ^i k ^m+k-i' > > and <S , ^^’”* 1 fc_1J ^n 
E+r l r iJ,) = 0 for 3 + 3', which implies 

f J’W-’li-l) = / , 

m+k — 1 m+k — 1 

for each fixed sequence ii, i 2 , • • •, 4-i- 

Let h = {n — 1) Sp(n), and 5' = {5/2) h , then (1211) holds. The proof is completed. ■ 

Now we come to the proof of Theorem [6l 
Proof of Theorem [6] 
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Denote y k = [x(k),x(k — !),••• , x(k — r)] T , then we have: 


V k - 1 C k y k) 


(27) 


where 


C k 


A°(k) A l (k) 

I 0 

0 I 


A T ~ l (k) A T (k) 
0 0 

0 0 


0 0 ••• / 0 

Let B k , = B k , and T' k , = Twh- Then {B' k . T' k } also forms an adapted sequence, 

and E {B' m+l \ T' m } has a 5-spanning tree. From Lemma QT) there exists 5' > 0 such that 

m+(n—1) Sp(n) 

IP | ^2 B'k has a ^'-spanning tree| T' m | > 5'. 

k=m -\-1 

Since 


diag(A°(A:)) 

0 •• 

• 0 

0 ’ 


diag(2L u (fc)) 

A 1 (k) • 

•• A T ~ 1 (k) 

1- 

ST 

F 

/ 

0 •• 

• 0 

0 


0 

0 

0 

0 

0 

I ■■ 

• 0 

0 

+ 

0 

0 

0 

0 

0 

0 •• 

■ I 

0 


0 

0 

0 

0 


> 5Mq + D k > 5(M 0 + D k ), 


where D k refers to the second matrix on the righthand side of the first line. From Lemma [7] 
Corollary IU and Remark [7] Remark [8] there exist 5"(p) 6 (0,5'), p > 1 such that 

m+p(n— 1) Sp(ra) 


P 

k=m -\-1 

■V _ T-rfc'(n-l)Sp(n) 


t r" — rr 

^k' 1 lfc=(fc' — l)(n—1) Sp(n)+1 

sequence and 


{ Iib-ryyu,— A-J x 

J] C k is 5"(p)-SIA,p =1,2,--- I T' m \ > 5'. 

L—<m_Ll ' 

Ck , T 7 " = J^, (n _ 1)Sp(n) , then forms an adapted 


V 

p { n is ^"(Pt-SIA ,P = 1,2, • • ■ I K] > S'. 


k =1 
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Let k n = n(n) + 1, then from Lemma [4] and Lemma [9], 

kn 

P { W C 'm+i is S " (I)' fen -scrambling| J'"} 

i —1 
V 

> P { II C 'm'+k is 6"(p)~ SIA, m' = m, ■ ■ ■ ,m + k n — l,p = 1,2, ■ ■ ■ ,| 

k k =1 

> <5 ,a A 

Let C" = ri/”=( m -i)fc„+i and K, = K,k„ - Th en {C", Ff) forms an adapted sequence, and 

P{C'" +1 is 5" (1) fcn - scrambling | F’f] > 5 ,k \ 

From the Second Borel-Cantelli Lemma (Lemma [8]), this implies 

P{C" is -scrambling infinitely often} = 1. 

From Lemma [3] and Remark HI we have 

p {^ A (n^)=o}=i. 

m= 1 

This implies the system (ITT!) reaches consensus almost surely. On the other hand it is not difficult 
to show that consensus of (1271) can imply consensus of (l20l) . And the proof is completed. 

From Theorem [6} we can have the following two corollaries. The first one is for deterministic 
case. 

Corollary 2: Let {A°(k) : ■ ■ ■ , A T (k)} be a deterministic sequence with A°(k) > 51 for some 
5 > 0, and [A l {k)\jj = 0 for i = 1,2, ■ ■ ■ ,r, j = 1 ,n. If there exists h > 0 such that 

Ylk=m+i -S(&) has a 5-spanning tree, then the system (l20l) will reach consensus. 

The second one is for a random sequence that will become independent after a fixed length 
of time. 

Corollary 3: Let {A°(k),--- ,A T (k)} be a random sequence with A°(k) > 51 for some 
5 > 0, and [A l (k)]jj — 0 for i — 1, 2, • • • , r, j — 1, • • • , n. And there exists > 0 such that 
A°(k), ■ ■ ■ , A T {k ) is independent of A°(k’), • • • , A T (k') whenever k! > k + N. If there exists 
h > 0 such that YHk=m +i p B{k) has a 5-spanning tree, then the system (l20l) will reach consensus 
almost surely. 

Proof: Let Fu = a{A°(m), • • • , A r (m), m < k} be the cr-algebra formed by A°(m), • • • , A T (m), 
m < k. Then {A°(k),--- , A T (k), Fk} forms an adapted process. From the assumption on 
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{A°(k),--- ,A T (k)}, A°(k'),--- ,A T (k') is independent of Bp whenever k! > k + N. Since 
B(k ) > 0, we have 

k-\-N-\-h k+N+h k+N+h k+N+h 

E{ Y B{m)\B k }>E{ Y B{m)\B k } = E{ Y 5 M} = 5^ EB{m). 

m=k-\- 1 m=k-\-N-\-l m=k-\-N-\-l m=k-\-N-\-l 

Thus, Birn^Bk} has a 5-spanning tree and the conclusion follows from Theorem 

m ■ 
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